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Abstract 

We conducted wide-field JHK S imaging polarimetry toward NGC 2024, which is a massive star-forming 
region in the Orion B cloud. We found a prominent and extended polarized nebula over NGC 2024, and 
constrained the location of illuminating source of the nebula through the analysis of polarization vectors. 
A massive star, IRS 2b with a spectral type of 08 — B2, is located in the center of the symmetric vector 
pattern. Five small polarized nebulae associated with YSOs are discovered on our polarization images. 
These nebulae are responsible for the structures of circumstellar matter (i.e., disk/envelope systems) that 
produce strongly polarized light through dust scattering. For the point-like sources, we performed software 
aperture polarimetry in order to measure integrated polarizations, with which we detected candidate 
sources associated with circumstellar material by selecting sources with a larger polarization than that 
estimated from the extinction of foreground material. We investigated the fraction of highly polarized 
sources against the intrinsic luminosity of stars (cx mass) , and found that the source detection rate remains 
constant from low (brown dwarfs) to higher luminosity (solar- type) stars. This result indicates the relative 
disk scale-hight is rather independent of the stellar mass. We confirmed the result using our polarimetry of 
the stars with known spectral types in NGC 2024. We found five young brown dwarfs with highly polarized 
integrated emission. These sources serve as direct evidence for the existence of disk/envelope system around 
brown dwarfs. We investigated the magnetic field structure of NGC 2024 through the measurements of 
dichroic polarization. The average position angle of projected magnetic fields across the region is found to 
be 110°. We found a good consistency in magnetic field structures obtained using near-infrared dichroic 
polarization and sub-mm/far-infrarcd dust emission polarization, indicating that the dichroic polarizations 
at near-infrared wavelengths trace magnetic field structures inside dense (Ay < 50 mag) molecular clouds. 

Key words: circumstellar matter — infrared: stars — ISM: individual (NGC 2024) — polarization — 
stars: formation 



1. INTRODUCTION 

NGC 2024 is a massive star forming region in the Orion B giant molecular cloud at a distance of 415 pc (Anthony- 
Twarog 1982). In the optical images, NGC 2024 shows a bright nebulosity with a prominent central dark dust lane 
located in front of the ionized gas. Imaging observations at the near-infrared wavelengths clearly revealed an embedded 
cluster in this region (e.g., Barnes et al. 1989; Lada et al. 1991; Haisch ct al. 2001), which is high in stellar density 
(~ 400 stars pc -3 , Lada 1999) and is young in mean age (0.3 Myr, Meyer 1996). The structure of the HII region 
toward NGC 2024, forming at the edge of a molecular cloud, was studied in detail (Crutcher et al. 1986; Barnes et 
al. 1989). Though the dominant source of ionization powering the HII region has not been identified until recently 
due to heavy obscuration by dust around the central region, Bik et al. (2003) suggested that the dominant source is 
IRS 2b located 5" north-west of IRS 2. They found that the spectral type of IRS 2b is in the range 08 V - B2 V, 
which is consistent with the intensities of radio continuum emission and recombination lines observed toward the HII 
region (Krugcl et al. 1982; Barnes ct al. 1989). NGC 2024 includes a number of protostars forming along the "star 
forming ridge" , a filamentary shaped dense molecular material, extending in the north-south direction. In the ridge, 



2 



Kandori et al. 



dense dust condensations, FIR 1-7, were reported based on the 1.3 mm and 870 /mi observations (Mezger et al. 1988, 
1992). 

The magnetic field structure of NGC 2024 was studied by observing linearly polarized thermal emission from aligned 
dust grains (B± measurements) and the Zeeman splitting of molecular or atomic lines (B\\ measurements). Hildcbrand 
et al. (1995) and Dotson et al. (2000) reported far-infrared (100 /mi) polarization map covering ~ 3' x 3' toward NGC 
2024. They found that the overall magnetic field orientation is in the east-west direction. There is a local magnetic field 
structure at and around the star forming ridge, but its detailed structure could not be resolved on their polarization 
vector map. Matthews, Fiege, & Moriarty-Schieven (2002) studied the magnetic field structure of the star forming 
ridge in NGC 2024 based on the sub-mm dust continuum observations using the polarimetry mode of JCMT/SCUBA. 
They revealed the detailed magnetic field structure along the ridge, which is consistent with the previous 100 /tm 
lower resolution polarization map. On the basis of their polarization map at 850 /tm, they modeled the magnetic field 
structure of the star forming ridge as helical magnetic fields surrounding a curved filamentary cloud, and as magnetic 
field swept by the ionization front of the expanding HII region. 

Crutcher et al. (1999) performed circular polarimetry toward NGC 2024 with the Very Large Array (VLA) to 
measure the Zeeman splitting of the absorption lines of HI and OH. They found Bn smoothly varies from fiG in the 
northeast to ~ 100/iG in the southwest of the star forming ridge, and suggested that the orientation of magnetic fields 
toward the lines of sight varies across the region. 

At the near-infrared wavelengths, polarimetry of a few bright stars in NGC 2024 were reported (IRS 1: Wilking et 
al. 1980, IRS 1,2: Heckert &; Zeilik 1981). More recently, Moore & Yamashita (1994) conducted imaging polarimetry 
toward the infrared nebula associated with FIR 4. However, the field of view of their map was too small (~ 40" x 40") 
to reveal polarization distribution over the entire surface of NGC 2024. 

In this paper, we present the polarization images of the NGC 2024 star forming region, which is a part of our ongoing 
project of JHK S polarimetry of star forming regions. Our observations are sufficiently deep (J = 19.2 mag, S/N = 5) 
and wide (7.7' x 7.7') covering a large extent of NGC 2024. Near-infrared imaging polarimetry is particularly useful 
to trace polarized light from star forming regions compared with that at optical wavelengths where dust extinction is 
severe. Our wide-field JHK S polarization images can be used to reveal large-scale infrared reflection nebulae (IRNc); 
to detect smaller-scale IRNe associated with young stars which is responsible for the circumstellar material (e.g., 
disk/envelope system), producing strongly polarized light through dust scattering; and to delineate the magnetic field 
structure through measurements of dichroic polarization of point sources. 

In §2 we describe the observations and data reduction. We derived the polarization maps of NGC 2024 and 
constrained the illuminating source of the large-scale IRNe through the polarization vector analysis in §3.1. In §3.2 
we present the newly found small-scale IRNe associated with young stars, and discuss their physical properties. In 
§3.3 we present the results of software aperture polarimetry of point sources. The polarization properties of point 
sources are investigated based on the measurements of both polarization degrees and colors due to dust extinction. In 
addition, we reveal the magnetic field structure over NGC 2024 obtained from the polarization angle of each source, 
and compare it with the data taken at other wavelengths. Our conclusions are summarized in §4. 

2. OBSERVATIONS AND DATA REDUCTION 

We carried out near-infrared polarimetric observations of NGC 2024 using the imaging polarimeter SIRPOL (po- 
larimetry mode of the SIRIUS camera: Kandori et al. 2006) on the 1.4 m telescope IRSF (Infrared Survey Facility) at 
the South African Astronomical Observatory (SAAO). SIRPOL consists of a single-beam polarimeter (rotating achro- 
matic [1 — 2.5/un] wave plate and high-extinction-ratio wiregrid analyzer) and the Jffif s -simultaneous imaging camera 
SIRIUS (Nagayama 2003), which has three 1024 x 1024 HgCdTe infrared detectors (HAWAII array). IRSF/SIRPOL 
enables deep and wide-field (7.7' x 7.7' with a scale of 0.45" pixel -1 ) imaging polarimetry at JHK S simultaneously. 
Polarizations of all the 2MASS-detected sources (15.8, 15.1, and 14.1 mag at J, H, and K s , S/N = 10) can be measured 
within the accuracy of 1 %. 

Observations were made on the night of 2006 February 13. The 1.25 /tm (J), 1.63 /mi (iJ), and 2.14 /mi (K s ) 
imaging polarimetry data of NGC 2024 were obtained simultaneously. The 10 sec exposures at 4 wave plate angles (in 
the sequence of 0°, 45°, 22.5°, and 67.5°) were performed at 10 dithered positions (1 set). We repeated the same set 
9 times toward the object. The sky frame were taken in the same manner. The total integration time toward NGC 
2024 was 900 seconds per wave plate angle. The typical seeing during the observations was ~ 1.4" (3 pixels) at J. We 
obtained twilight flat at the beginning and end of nights. In order to make good flat frames, we combined 23 sets of 
twilight flat observations (17 nights) taken in the period from December 2005 to February 2006. 

We processed the observed data in the same manner as described in Kandori et al. (2006) using the IRAF 1 (flat- 
field correction, median sky subtraction, and frame registration). The Stokes /, Q, and U images are obtained using 
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1 = (Ia° +^45° +^22.5° +^67.5°)/2, Q = Iq° — 1^5° , and U = 1-22.5° ~ -^67.5° ■ The image of the polarization degree P 
and polarization angle 9 can be derived by P = \/ Q 2 + U 2 /I and = 0.5 x atan(f//Q). We calibrated I sky level to 
be consistent with the surface brightness derived using 2MASS images covering a much larger field-of-view. The 10cr 
limiting magnitudes for surface brightness of I were 19.6, 19.3, and 17.3 mag arcsec -2 for J, H, and K s , respectively 
In Fig.l we show images of the JiJET s -composite intensity (I) and polarized intensity (PI). We note that the PSF 
of unpolarizcd stars are not perfectly canceled on the PI image, because the seeing size changes slightly during the 
observations to take Q and U. 

For the source detection and photometry on Stokes I images, we used the IRAF daophot package (Stetson 1987). 
We detected stars having a peak intensity greater than 10cr above local sky background. The overlooked sources in 
automatic detection were collected by eye inspection. The detected number of sources was 240, 365, and 376 at J, 
H, and K s , respectively We measured instrumental magnitudes of stars with aperture photometry. The aperture 
radius was 3 pixels. The number of stars detected in all three bands with a photometric error of less than 0.1 mag is 
211. The lOer limiting magnitudes were 18.5, 18.0, and 16.5 mag for J, H, and K s , respectively. In order to obtain 
plate solutions, we matched pixel coordinates of detected sources and celestial coordinates of their counterparts in the 
2MASS Point-Source Catalog (Cutri et al. 2003), and applied the IRAF imcoords package to the matched list. The 
coordinate transformation error is ~ 0.1" (rms). We performed photometric calibration using the 2MASS-PSC. The 
magnitude and color of our photometry were transformed into the 2MASS system using, 

MAG 2 mass = MAGirsf + a t x COLOR msF + (1) 
COLOR 2 MAss = «2 x COLOR IRSF + p 3 , (2) 

where the coefficients a.\ are 0.029, —0.024, and 0.008 for J, H, and K s , respectively. The coefficients 0,2 for J — H 
and H — K s are 1.032 and 0.997, respectively. We note that the coefficients (3i and /3a include both zero point and 
aperture correction. 

We carried out software aperture polarimetry of point sources on the combined intensity images for each wave 
plate angle (Io°, 122.5°, 1^5°, and 167.5° )■ This is because the center of point sources (i.e., aperture center) cannot be 
determined satisfactory on the Q and U images. We discuss the results in the H band where nebula contamination is 
less than in the J band and scattering efficiency is higher than in the K s band. Based on aperture photometry data of 
each angle image, the Stokes parameters for each source were derived as described above. Point sources on each angle 
image have a slightly different PSF due to varying seeing size (~ 1.3" — 1.5" at J). We thus used different apertures 
for each image in order to measure the same fraction of stellar flux falling in the aperture. The adopted aperture 
sizes were 3.00, 3.24, 3.08, and 3.27 pixels for Iqo, ^22.5°, -^45°, and i67.5°; respectively. The sky annulus was set to 10 
pixels with 5 pixel width. Since polarization degree P is a positive quantity, derived P values tend to be overestimated 

especially for low S/N sources. To correct the bias we calculated the debiased P using P db = VP 2 -SP 2 , where SP 
is the error in P (Wardle & Kronberg 1974). 

3. RESULT AND DISCUSSION 

3.1. large- Scale Infrared Reflection Nebula (IRN) 

Although near-infrared (H and K) polarimetry around a small nebula associated with FIR 4 was reported previ- 
ously (Moore & Yamashita 1995), no wide-field polarimetry of NGC 2024 has been carried out. Our wide and deep 
observations with SIRPOL revealed the distribution of near-infrared polarization in NGC 2024 for the first time. In 
Fig. 2, we show the /, PI, P, and 8 of NGC 2024 in the H band. The central dark lane running north-south is located 
foreground of the nebula and obscures the light from H II region (Fig. 2: upper-left). In the same panel, we show 
the location of far-infrared sources (FIR 1-7, Mczger et al. 1988, 1992) associated with a filament of dense material, 
"star-forming ridge", located just behind the H II region. Among the far-infrared sources, FIR 4 is associated with a 

2 /im source (Moore & Chandler 1989) as well as reflection nebula (Moore & Yamashita 1995), both of which are also 
confirmed with our data (see §3.2). In addition to the large scale IRN, there are smaller scale IRN associated with 
stars, which we discuss in the next section. We found prominent and extended •polarization nebulosities over NGC 
2024 for the first time (Fig. 2: upper-right). In contrast to the other regions, e.g., M 42, the exciting source of NGC 
2024 is not directly identified due to the heavy obscuration by the central dust lane. Although IRS 2 was thought to 
be the exciting source of NGC 2024 until recently (e.g., Grasdalen 1974), Bik et al. (2003) proposed the the dominant 
ionizing source of NGC 2024 is IRS 2b located 5" north-west of IRS 2 (Fig. 2: upper-left). 

Our imaging polarimetry data can be used to constrain the location of ionizing sourcc(s) of the nebula through the 
polarization vector analysis assuming that the ionizing source and the illuminating source are the same. The vector 
map is shown in Fig. 3. We found the vector patterns clearly appear centrosymmetric, and normal lines of each 
vector intersect at the position near the center of NGC 2024 shown by white plus symbol. The location of dominant 
illuminating source(s) of NGC 2024 is then determined to be located within the white circle (la error circle) in Fig. 3. 
We note that the locations of both IRS 2 and IRS 2b arc included in the error circle but IRS 2b is closer to the center. 
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The vector patterns at J, H , and K s are consistent with each other and result in the coincident symmetry center 
position. We colored polarization vectors black or green based on the angle 9 between each vector and the line toward 
white plus symbol; vectors which seem to follow centrosymmetry, showing 9 = 90 ± 20°, arc in white, and other vectors 
are in green. Polarization vectors are mostly colored in white, indicating that most of the large-scale nebulosities arc 
illuminated from star(s) in the central region enclosed by white circle. The green vectors tend to be distributed over 
the north-south dark dust lane and around small-scale IRNe, which can be explained by superposition of the large-scale 
nebulosities and extra polarization fields, i.e., dichroic absorption of foreground materials and intrinsic polarization of 
small-scale IRNe. 

3.2. Small- Scale IRN 

We found five new small-scale IRNe (IRN 1-5) associated with stars that are shown in the PI image in Fig. 2. 
An object enclosed by dashed line circle is the IRN associated with FIR 4, which was studied previously by Moore 
& Yamashita (1995) with near-infrared polarimctry. We show PI and vector maps of IRN in Fig. 4. It is clear 
that the pattern of the polarization vector appears centrosymmetric indicating that the nebulosity is associated with 
the central star, and responsible for the structure of circumstcllar matter (e.g., disk/envelope systems) that produces 
strongly polarized light through dust scattering. The small nebulae in a larger scale nebulosity or embedded in the 
cloud should be treated with caution, because the polarization toward such nebulae is a superposition of different 
polarization vectors (i.e., scattered light polarization and dichroic polarization). We subtracted local background 
(constant value) around small-scale IRNe from the Stokes Q and U images in order to suppress the contribution from 
larger nebulosities. We determined the Q and U background levels at the position 20" away from the center of IRN, 
and simply subtracted it from the original. We did not take the effect of dichroic polarization into account, because 
the resulting vector pattern of IRNe after the subtraction of local background was reasonably centrosymmetric and 
the estimate of dichroic polarization is not straightforward. 

All IRNe have certain structures on the PI images. The spatial extent of small-scale IRNe is typically about 
5000 — 10000 AU on the PI images, comparable to the size of "envelopes" rather "disks" . IRNe show an elongated 
polarized emission pattern whose position angle is roughly 90, 135, 30, 45, and 135° for IRN 1-5, respectively. The 
emission pattern of IRN 1, 2, and 4 is clearly butterfly-shaped, and the dark lane extending perpendicular to the 
bright region can be seen in these IRNe. For IRN 3 and 5, the elongated emission pattern can only be seen on one side 
of the central star. These characteristics agree with the polarization picture of disk/envelope system around young 
stars (e.g., Nagata, Sato, & Kobayashi 1983; Hodapp 1984; Sato et al. 1985). The radiation from the star causes 
heavy extinction toward the direction of disk (seen as a dark lane), but can escape toward the polar directions and 
then scattered by dust at the surface of the envelope or the cavity created by the outflow (seen as a butterfly-shaped 
or monopolar IRN). 

In Fig. 5, we show PI and vector maps of the IRN associated with FIR 4 in the J, H, and K s bands. The 
counterpart of FIR 4 is best seen in the K s band and located at (AR.A., ADcc.) = (7.2", —7.2") in Fig. 5. Our vector 
maps show good agreement with those in the H and K bands reported by Moore & Yamashita (1995). The IRN 
appears as a monopolar shape with position angle of ~ 135°, and shows a centrosymmetric vector pattern centered on 
FIR 4. The axis of IRN is consistent with that of the molecular outflow associated with FIR 4 (Chandler & Carlstrom 
1996). The position angle of the polarization vector toward FIR 4 is ~ 45°, perpendicular to the axis of IRN, which 
is the general tendency in young stellar objects with IRN (see above). 

In Tabic 1, we summarize the observational properties of IRNe and its central stars. The polarization degree of each 
IRN (column 4-6th) is not large, because unpolarizcd light from the HII region (e.g., free- free ionized gas emission) 
dilute the polarized light. If we subtract the background intensity around IRN, the polarization degrees increase 
to the level comparable with typical reflection nebulae in other star forming regions as shown in parentheses. The 
IRN 1-5 are associated with stars IRS 12, 4, 3, 28 and 30, respectively. From their near-infrared colors, we found 
that the central stars of five IRNe are young, especially in a pre-main-sequence phase. The infrared color of the five 
sources were investigated in the JHKLN bands by Haisch et al. (2001). Based on the 2.2— 10 (im spectral index 
a (= d\og\F x /d\og\) measurements, the five sources, IRS 12 (IRN 1), IRS 4 (IRN 2), IRS 3 (IRN 3), IRS28 (IRN 
4), and IRS 30 (IRN 5), were classified into Class II, Class I, Class II, Class II, and Flat Spectrum, respectively. 
This is consistent with the fact that they are accompanied by nebulous circumstcllar materials. Following the method 
described in §§ 2, we performed aperture polarimetry of the central stars (column 12-13). Except for the saturated 
source IRS 12 and the source IRS 3 with faint IRN, we found that the polarization vector angle 9 of stars is roughly 
perpendicular to the major axis of IRN (column 3). The characteristic is in good agreement with those for young 
stellar objects accompanied by IRNe. These sources serve as good target for future studies in higher resolution at 
various wavelengths (e.g., large Opt-IR telescopes, ALMA) as well as in theoretical modeling. 

3.3. Aperture Polarimetry of Point-like Sources 

We performed software aperture polarimetry of point-like sources in order to measure the integrated polarization of 
possible unresolved nebulae. Our photometry and polarimetry of 211 stars with JHK S detections are shown in Table 
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2. The polarization vectors in the J, H, and K s bands are quite consistent with each other. For the sources with 50\ 
of less than 3° (63 sources), the correlation coefficient for (Oh, @j) and (Oh, 0k b ) is 0.99 and 0.98, respectively. The 
mean of the difference in 0\ is small, which is —1.4° for < Oh — 0j > and —1.7° for < Oh — 0j< s >. For the sources 
with P\/8P\ of greater than 10 (64 sources), the linearly fitted slope of Ph vs. Pj and Ph vs. Pk s diagram is 1.52 
and 0.50, with the correlation coefficient of 0.98 and 0.91, respectively. These slopes are consistent with the empirical 
value of 1.61 and 0.61 from the relation P oc A _Q , where a = 1.8 ±0.2 (Whittet 1992). In the following, we discuss 
the aperture polarimetry results in the H band where nebula contamination is less than in the J band and scattering 
efficiency is higher than in the K s band. 

3.3.1. Highly Polarized Sources (HPS) 

We can detect candidate stars with intrinsic polarization (i.e., with circumstellar material) by selecting the sources 
with larger polarizations than values estimated from the dichroic absorption due to foreground material. Firstly, we 
made two color diagram for all the sources detected in the JHK S bands in order to classify the sample into three 
groups of the dwarf+giant, PMS stars, and protostars. The result is shown in Fig. 6a. We found 143 stars in the 
dwarf+giant star region (filled circle), 63 stars in the PMS star region (open diamond), and 5 stars located on right 
to the PMS region, i.e., protostar region (open circle). We plotted these sources in H — K s — 0.15 versus Ph diagram 
shown in Fig. 6b. The value of 0.15 mag is the mean intrinsic color of stars (mostly G to M dwarfs) toward the 
object calculated using the Galaxy model (Wainscoat & Cohen, 1992; Cohen, 1994), so that the horizontal axis of 
Fig. 6b corresponds to the mean Eh-k b - The dashed line denotes the mean relationship between extinction and 
polarization for interstellar molecular clouds, Pk(%) — 2.23r^/ 4 , from Jones (1989). A linear least absolute deviation 
fitting to our H — K s — 0.15 and Ph data resulted in Ph = 4.21 x Eh-k s = 1-69 x Ah (solid line), where we forced 
the intercept value to be zero in the fitting. Our fitting result shows a similar relationship to the relation from Jones 
(1989). The polarization of stars with Ph/Ah ~ 1-69 is dominated by dichroic polarization produced by aligned 
interstellar dust grains. Observational upper limit of the relationship (dot-dashed line) can be described using the 
equation P max = tanh(rK ) where tk = (1 — v) t k / (1 + v) when the parameter r\ is set to 0.875 (Jones 1989). 

The sources with larger Ph/Ah values can be better candidates to show excess (intrinsic) polarizations, and hense 
good targets for further observations in order to study circumstellar materials around stars. We calculated Ph/Ah 
for each star, where Ah was measured on the two color diagram (Fig. 6a). The histogram of Ph/Ah is shown 
in Fig. 7. In the histogram, there is a clear peak around Ph/Ah ~ 1.69, which corresponds to the sources mainly 
affected by interstellar polarization. We found that the dwarf+giant stars (filled circle) and PMS stars (open diamond) 
show similar distribution in Fig. 6b. If we regard the stars with Ph/Ah > 3.38, two times larger Ph/Ah than the 
interstellar polarization, as highly polarized sources (HPS), the fraction of HPS is 32 %, 31 %, and 35 % for all the 
stars, dwarf+giant stars, and, PMS stars, respectively The protostar candidates are not included in this statistics, 
because their intrinsic colors are not known. It is an interesting result that there is no tendency of higher polarization 
in earlier evolutionary stages or larger infrared excess. However, since the part of the locus of classical T Tauri stars 
overlaps with the reddening band, some fraction of PMS stars can be included in the sample of dwarf+giant stars. 
Thus the HPS fraction for the dwarf+giant group may cause significant contamination from PMS stars and caution 
must be applied. 

3.3.2. Stellar Mass versus P H /A H 

It is of particular interest to investigate the fraction of HPS against the luminosity of stars (oc mass), because the 
information should be closely related to the dependence of the population of stars with circumstellar disk on the mass 
of stars. Though the probability of disk inclination, orientation, and suffered extinction from foreground materials is 
expected to be random, other factors, the shape (scale height) and evolutional timcscalc (lifetime) of disk, may have 
some dependence on stellar mass. 

It is likely that the effect of disk lifetime on the detection ratio of HPS is not significant. As shown in Fig. 6a, many 
sources are located in the Class II region in the two color diagram. By using J — H vs. K — L diagram, Haisch et al. 
(2000) reported that infrared excess stars (Class II and protostars) occupy ~ 80 % fraction of the stellar population in 
this region. Since the age of most sources could be younger than or comparable with disk lifetime (~ 10 Myr, Strom 
& Edwards 1993), the observed HPS fraction should not be affected by the disk lifetime. 

Theoretically very low mass stars (i.e., brown dwarfs) are expected to have a more flared (larger scale height) 
disk than those for more massive stars due to its small gravity (Walker et al. 2004). Highly flared disks result in a 
large fraction of obscured stars due to extinction, so that only scattered light is observed, making the source a HPS. 
Considering the disk optical depth (i.e., size and mass) and its flaring for several models, the disk obscuring probability 
of brown dwarfs in the K band (2.2/mi) was estimated to be larger than that for Class II stars by a factor of ~ 1 — 3 
(Walker et al. 2004). Thus, a stellar mass dependence of HPS detection rate may be observed. 

To derive intrinsic luminosities, we performed dereddening of stars by subtracting Ah from observed H magnitudes. 
We then derived the fraction of sources with Ph/Ah > 3.38 against dereddened H mag at each bin separated by 1 mag. 
The result is shown in Table 3. Note that stars with H > 13.3 mag are brown dwarfs at the employed distance of 415 
pc (Anthony- Twarog 1982). We found that the distribution of polarization fraction remains constant from low (brown 
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dwarfs) to higher luminosity (solar-type) stars, although the error in polarization fraction is not small. According to 
the argument above, our finding suggests that the dependence of disk scale height on the mass of central star is not 
significant in spite of the theoretical model that predicts more flared disks for lower mass stars. It is noteworthy that 
a very similar result is obtained in the OMC-1 region (Kusakabc ct al. 2006). Quantitative studies on the origin of 
flat HPS distributions will be discussed in a later paper. 

Recently, Levine et al. (2004) reported the spectroscopic studies of low-mass stars and brown dwarfs distributed 
over 10' x 10' region in NGC 2024. Their list of 71 stars with known spectral type can be more useful to investigate 
the relationship between stellar mass and Ph /Ah ■ We matched the coordinates of their sources with our polarimetry, 
and found 49 pairs. We show the stellar mass versus Ph/Ah diagram in Fig. 8. The vertical and horizontal broken 
lines denote the upper limit of brown dwarf mass (0.08 M@) and Ph/Ah = 3.38, respectively It is clear that there is 
no obvious tendency for higher Ph/Ah for lower mass stars. If we divide the sources into two groups at the line of 
Al = 0.08M Q , the median Ph/Ah value with (rmsj/^/n error is 2.11 ±0.53 and 1.62 ±0.43 for the lower and higher 
mass groups, respectively. There can be a tendency of slightly higher Ph/Ah for brown dwarfs than that for solar-type 
stars, but, within la, those two values are identical. We should also note the difference in suffered mean extinction 
for both groups is small (< Ah >brown dwarf= 1-5 mag and < Ah >soair-typo= 1-9 mag). We conclude that the stellar 
mass vs. Ph /Ah relationship is generally flat in the mass range from brown dwarfs to solar-type stars. 

In Fig. 8, the sources located above the horizontal broken line can be classified as HPS. Their properties are 
summarized in Table. 4 2 . We found ten HPSs which consist of five brown dwarfs and five low-mass stars. It is most 
likely that these sources are associated with circumstellar materials which produce highly polarized light through dust 
scattering. The existence of circumstellar disk around young brown dwarfs has been indirectly suggested (e.g., infrared 
excess emission, Muench et al. 2001). In addition to a few samples revealed by imaging of disks using HST (Bally 
ct al. 2000) and outflow detection with spectro-astrometry (Whelan et al. 2005), our five samples serve as "direct" 
evidence for the existence of disk/envelope system around brown dwarfs (see also, Tamura et al. 2006; Kusakabe et 
al. 2006). 

3.3.3. Magnetic Field Structure 

The aperture polarimetry of stars provides important information on the direction of magnetic fields. If we assume 
the normal grain alignment, i.e., the spin axis of elongated dust grain align parallel to the magnetic fields (Davis 
& Greenstein 1951), the direction of magnetic fields projected onto the sky can be inferred from the direction of 
polarization vectors of stars (B \\ E, e.g., Wcintraub, Goodman, & Akeson 2000). Though a grain alignment mechanism 
caused by gas streaming motion (Gold 1952) may be the dominant process near or around sources with molecular 
outflows, we believe our assumption is plausible, because we discuss the structure of magnetic fields not only around 
each outflow sources but also across the entire surface of NGC 2024. 

In Fig. 9, we show the histogram of polarization position angle, 9h, of stars in NGC 2024. It is clear that the stars 
with Ph > 3 % show a peak at the position angle of 110°. The large scale distribution of magnetic field in NGC 2024 
is thus expected to be in the direction of 110°. In the right panel, the stars with low polarization (Ph < 3 %) shows 
less clear tendency. In the histogram, we show the number of protostars, PMS stars, and dwarf±giant stars in dark 
gray, gray and white colors, respectively. There is no large difference in the shape of histograms among three groups 
with different evolutionary stage. 

We show the distribution of polarization vectors for each star in Fig. 10. Though the overall distribution of 
polarization angle is ~ 110°, we found a bend of magnetic fields around the northwest of the star-forming ridge of 
NGC 2024. Since the bend of magnetic field orientation is localized on the polarization map, it is most likely that 
the structure is not associated with ambient material around NGC 2024 but associated with a smaller-scale structure 
inside the region, i.e., the star forming ridge. 

Matthews, Fiege, & Moriarty-Schieven (2002) studied the magnetic field structure of NGC 2024 based on the sub- 
mm dust continuum observations using the polarimetry mode of JCMT/SCUBA. We overlaid the polarization vector 
of stars on their polarization map in Fig. 11. Vectors from the sub- mm emission polarimetry (blue lines) were rotated 
by 90 degrees so as to be compared with the E- vector of dichroic polarization in the H band (yellow lines) . We hither 
plotted rotated (90 degrees) polarization vectors taken with 100 fim dust emission (red lines, Hildebrand et al. 1995; 
Dotson et al. 2000). The background image in Fig. 11 is the 850 /jm dust continuum intensity map taken with 
JCMT/SCUBA. The gray broken line shows the contour of 0.9 Jy beam -1 , which corresponds to Ay ~ 50 mag when 
we use the dust opacity k 85 o ,,„, = 0.02 cm 2 g _1 and dust temperature Td = 18 K. The region enclosed by the gray 
broken line is too highly obscured to be traced at the near-infrared wavelengths. 

We found that the direction of magnetic field obtained from both dust emission (blue and red lines) and extinction 
(yellow lines) is generally consistent with each other in the region around and outside the contour of Ay ~ 50 mag. The 
agreement between dichroic and dust emission polarizations suggests that dichroic polarization traces well the magnetic 
field structure if lying dense dusty material is less than or as much as its penetration depth. On the other hand, dichroic 
polarization vectors do not show good correlation with those obtained from dust emission inside the contour of Ay ~ 50 



Full listing of all sources is shown in Table 5 (electric version only). 
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mag. It is a natural result if we consider that the dichroic polarization in the H band only trace the magnetic field 
foreground to the star-forming ridge. The disagreement between dichroic and dust emission polarizations suggests 
that magnetic field structure toward the star-forming ridge traced using the dust emission polarimetry is localized 
in both the plane-of-sky and the line-of-sight direction. We note that the line-of-sight variation in magnetic field 
strength (B\\) toward the star-forming ridge is also observed by Crutcher et al. (1999) based on the measurements 
of the Zeeman splitting of the absorption lines of HI and OH with Very Large Array (VLA). On the basis of the 
sub-mm dust emission polarimetry data, Matthews, Fiege, & Moriarty-Schieven (2002) explained the magnetic field 
structure toward the star-forming ridge using two different models; one is helical magnetic fields surrounding a curved 
filamentary cloud, and another is magnetic fields swept by the ionization front of the expanding HII region. It is hard 
to differentiate the two models based on our near-infrared dichroic polarization data, because the data can not trace 
the dense region (Ay > 50 mag) of the star-forming ridge. 

It is noteworthy that our result is contrary to the view that the polarization of background starlight in the optical 
to near-infrared does not trace magnetic fields (i.e., grain alignment) inside cold and dense molecular clouds (e.g., 
Goodman et al. 1995). For example, Arce et al. (1998) suggested that the dichroic polarization of background stars 
toward the Taurus dark cloud only trace the region up to Ay ~ 1.3 mag, indicating that the dichroic polarization 
does not trace magnetic field inside cold dark clouds due to significant depolarization effect. Their criticism may not 
be the case for other star forming regions where the environment is different (e.g., higher temperature 3 ) from those 
in cold dark clouds. In addition to the NGC 2024 region, a good consistency in magnetic field structure obtained 
using dichroic polarization and dust continuum polarization is confirmed in the OMC-1 (Houde et al. 2004; Kusakabe 
et al. 2006) and NGC 2071 (Matthews, Fiege, & Moriarty-Schieven 2002; Tamura et al. 2006) region. We thus 
conclude that the dichroic polarization at near-infrared wavelengths is a good tracer of magnetic field structures than 
previously thought and wide-field near-infrared polarimetry is useful to reveal the magnetic field structures in various 
star-forming regioins. 

4. SUMMARY 

We conducted deep and wide-field JHK S imaging polarimetry toward NGC 2024, a massive star-forming region in 
the Orion B cloud. This is the first imaging polarimetry covering almost all of the HII region as well as embedded 
young cluster in NGC 2024. Our conclusions are summarized as follows. 

1. We found a prominent and extended infrared reflection nebula (IRN) over NGC 2024 on our polarization images. 
We constrained the location of illuminating source of the nebula through the analysis of polarization vectors. A 
massive star, IRS 2b with a spectral type of 08 — B2, is located in the center of the symmetric vector pattern. 

2. We discovered five small-scale IRNe associated with YSOs in our polarization images. The pattern of polarization 
vector around the stars appears ccntrosymmetric indicating that these IRNe are illuminated by each central star 
and are responsible for the structure of circumstellar material that produces strongly polarized light through 
dust scattering. The polarized intensity of each IRN shows a butterfly-shaped emission pattern extending over 
5000— 10000 AU, which agrees well with the picture of disk/envelope system around young stars. 

3. We performed software aperture polarimetry of the 211 point-like sources detected on the Stokes I images at 
JHK S . We cataloged 64 highly polarized sources (HPS) as candidates associated with circumstellar material by 
selecting the sources with stronger polarizations than estimated from dichroic extinction. 

4. We investigated the fraction of highly polarized sources against the intrinsic luminosity of stars (oc mass), and 
found that the source detection rate remains constant from low (brown dwarfs) to higher luminosity (solar- 
type) stars. Since most stars in the region are young enough to have circumstellar disks, possible differences 
in disk lifetime should not affect our statistics. Though a stellar-mass-dependence of disk scale height (flaring) 
is a possible factor to naturally change the detection rate of highly polarized source, our result indicates that 
the effect is not significant in spite of the theoretical prediction of more flared disks for lower mass stars. We 
confirmed the result using our polarimetry of the stars with known spectral types toward NGC 2024.. 

5. By comparing our polarimetry with the spectro-photometry by Levine et al. (2004), we found five brown dwarfs 
with highly polarized integrated emission. These sources serve as direct evidence of circumstellar material around 
brown dwarfs. 

6. We investigated the magnetic field structure of NGC 2024 through the measurements of dichroic polarization of 
point sources. The position angle of projected magnetic field across the region is generally found to be 110°. We 
also found a bending of magnetic fields around the northwest of the star-forming ridge of NGC 2024. 

7. In the region of Ay < 50 mag, we found good consistency in magnetic field structures obtained using dichroic 
polarization at near-infrared and dust continuum polarization at sub-mm/far-infrarcd wavelengths. Similar 

3 The temperature of gas and dust was measured as T^n < 20 K (Gaume, Johnston, & Wilson 1992; Mauersberger et al. 1992) and 
X<j ~ 16 — 18 K (16 K toward the star-forming ridge: Mezger et al. 1988, 18 K toward the Orion B South region: Johnstone et al. 
2006) in NGC 2024. 



8 



Kandori et al. 



results are also found in the OMC-1 and NGC 2071 regions. These results strongly indicate that the dichroic 
polarization at near-infrared wavelengths traces magnetic field structures inside dense molecular clouds. 
8. In the region of Ay > 50 mag over the star-forming ridge, the dichroic polarization docs not show good correlation 
with those obtained from dust emission. It is a natural result if we consider that the dichroic polarization at 
near-infrared wavelengths can only trace the alignment of foreground dust. The disagreement between dichroic 
and dust emission polarizations suggests that magnetic field structure toward the star-forming ridge traced using 
the dust emission polarimetry is localized in both the plane-of-sky and the line-of-sight direction. 
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Tabic 1: Properties of small-scale IRNe 



Nebulae Central Stars 
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Size and position angle of IRNe on the PI images. The size of the nebula having intensity greater than 10cr above background was measured as the size of the 
IRN. The position angle denotes the orientation of polarized nebulosity of each IRN. 

Maximum polarization degree of IRNe was measured after the subtraction of local background on Q and U images. If we further remove the contribution from 

large-scale nebula by subtracting the local background intensity (I), the polarization degree of IRNe increase as shown in parentheses. 

Source number in Table 2. 

Star-ID based on Barnes et al. (1989). 

Magnitudes taken from 2MASS Point Source Catalog (Cutri et al. 2003). 

Evolutionary state of the central stars (YSOs) based on the 2.2 — 10 /im spectral index measurements (Haisch et al. 2001). 
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Table 2: Polarizations of Point- like Sources in NGC 2024 
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(mag) 



Ph/A h c 
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(°) 
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31 


5±12.8 
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20 


50±3.6 
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1 
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2 
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20 
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1 
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35 
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1 
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21 
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2 
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67 
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05 
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14 
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4 
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02 
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2 
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3 


57±0. 


07 
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1 
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7±11.8 


2 


10±0. 


15 


163±2.1 


3 
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03 
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16 
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o 
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01 
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1 


57±0. 


50 
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14 
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2 
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00 
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4 
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44 
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03 
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1 
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20 
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4 


47±0. 


47 
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1 


99±0. 


12 
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1 


92±0. 


20 


63±3.0 


2 


,05±1. 


02 


106±12.8 


1 


,57±0. 


18 


75±3.3 



1.24 
1.67 
2.71 
1.64 
0.33 
2.23 
0.90 
1.51 
2.48 
3.78 
3.27 
2.93 
2.38 
1.41 
1.02 
1.95 
1.70 
4.23 
3.69 
2.62 
1.49 
1.65 
1.02 
2.17 
1.37 
1.49 
2.25 
0.72 
1.84 
3.14 
3.68 
0.59 
2.06 
2.59 
0.96 
1.09 
1.43 
1.48 
1.38 
0.98 
1.90 
0.00 
1.50 
0.62 
1.03 
2.16 
2.11 
2.24 
0.72 
1.20 
1.90 
1.57 
1.23 
0.00 
1.30 
1.88 
1.45 
1.35 
1.79 
1.44 
2.72 
1.24 
2.25 
2.98 
1.97 
2.47 
1.17 
4.35 
1.03 
2.04 



Source ID of highly polarized sources with Ph /Ah > 3.38. 

Equatorial coordinates (J2000). 

Ah estimated on the two color diagram. 

Errors are estimated using SPh ■ 



Kandori et al. 



Tabic 2: continued. 



No. 


HPS ID 




R.A. 


) 


Dec. 

(° ' ") 




J 

(mag) 


H 

(mag) 


K B 
(mag) 


Ph 

(%) 


A H 
(mag) 


Ph/Ah 
{% mag -1 ) 


S H 

(°) 
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5 


41 


1 

4;) 


31 
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53 


32 


7 
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30 
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Table 2: continued. 
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Table. 3: Fractions of Highly Polarized Sources 



H a 


Fraction'' 


Number'' 


Total Number c 


(mag) 


(%) 






9.5 


31 ±8.9 % 


12 


39 


10.5 


29 ± 8.2 % 


12 


42 


11.5 


31 ±9.5 % 


11 


35 


12.5 


36 ±11.3 % 


10 


28 


13.5 


39 ± 13.0 % 


9 


23 



Dereddened H magnitudes. Each bin is separated by 1 mag. 

The fraction and the number of stars with Ph/Ah > 3.38 in each bin. 

Total number of stars in each bin. 
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Table. 4: Highly Polarized Sources with Known Spectral Type 



Star Name" 


Pu b 


A„ c 


Ph/Ah 


9 H d 


M subclass 13 


Mass-' 




(%) 


(mag) 


(% mag" 1 ) 


O 




(M ) 


FLMNLJ0541328-0151271 (57) s 


4.22±0.41 


0.94 


4.49±0.43 


147±2.8 


8.50 


0.02 


FLMN_J0541506-0158041 (25) 


3.45±0.09 


0.61 


5.65±0.14 


88±0.7 


7.75 


0.03 


FLMN_J0541456-0157332 (29) 
FLMN_J0541497-0156198 (55) 9 


2.14±0.10 


0.20 


10.61±0.52 


102±1.4 


8.00 


0.03 


19.88±3.89 


5.10 


3.90±0.76 


79±5.5 


7.00 


0.04 


FLMN.J0541537-0157503 (52) 


6.07±0.36 


1.58 


3.85±0.23 


108±1.7 


6.50 


0.06 


FLMN_J0541523-0157123 (46) 


8.20±0.30 


1.31 


6.25±0.22 


94±1.0 


5.75 


0.11 


FLMN_J0541534-0156482 (51) 


4.20±0.28 


1.02 


4.12±0.27 


96±1.9 


4.00 


0.23 


FLMN.J0541314-0154347 (30) 


7.43±0.75 


2.08 


3.58±0.36 


130±2.9 


3.25 


0.35 


FLMNJ0541357-0153483 (61) 


14.45±0.85 


1.17 


12.39±0.73 


112±1.7 


2.75 


0.37 


FLMN.J0541344-0154409 (02) 


4.40±0.03 


0.83 


5.31±0.04 


109±0.2 


1.75 


0.72 



: ' e * Taken from Table 3 of Levine et al. (2004). The numbers in parentheses of column 1 denote the source serial number 

in their catalog. Ah was calculated by multiplying originally cataloged Ay by a factor 5.71. 
M This work. 

9 Not cataloged in Table 2 but have available Ph measurements. 



NGC2024: Intensity (JHK S ) Polarized Intensity (JHK S ) 




R.A (J2D0O) R.A <J2DWt 



Figure 1 — Left: three-color composite of JHK S intensity (I) images toward NGC 2024. Right: three-color composite of JHK S polarized intensity (PI) 
images toward the same region. We note that there are bad pixel clusters on the J image around the upper-right corner and the middle of right side, both 
of which are masked by hatched line. The images are in logarithmic scale. 
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Figure 2 — Upper- left: Intensity (I) image of NGC 2024 in the H band (logarithmic scale). Plus symbols denote 
far-infrared sources (FIR 1-6: Mezger et al. 1988; FIR 7: Mezger et al. 1992). Some of bright stars are marked with 
cataloged IRS numbers in Bernes (1989). Upper-right: Polarized intensity (PI) image of the same region 
(logarithmic scale). Newly found small infrared nebulae are enclosed by circles. A nebula enclosed by broken line 
circle was previously studied with infrared polarimetry (Moore & Yamashita 1995). Lower-left: Polarization degree 
(P) image. Ph is ~ 30 % at maximum in the H band. The regions with low S/N (I/5I < 10) are masked. 
Lower-right: Polarization angle (9jj) image. 




Figure 3 — Polarization vector map on the intensity image at H. White plus symbol indicate expected location of 
the illuminating source of NGC 2024. White circle around plus symbol denotes la error circle (see text). The 
polarization vectors on the low intensity region (I/5I < 10) are not plotted. 
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Figure 4 — Left column: Polarization vector maps on the / images of small IRNe in NGC 2024. The polarization 
vectors on the low intensity region (IRN 1-3: I /SI < 10, IRN 4-5: I /SI < 5) are not plotted. Right column: 
Polarized intensity images. All the images arc in the H band and in logarithmic scale. 
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Figure 4 — 



continued. 
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Figure 5 — Left column: Polarization vector maps superposed on the / images of 40" x 40" region around FIR 4. 
The polarization vectors on the low intensity region (I /SI < 10) are not plotted. Right column: Polarized intensity 
images. The images are in logarithmic scale. 
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Figure 6 — Left panel: Two color diagram of the stars toward NGC 2024. Filled circles, open diamonds, and open circles denote dwarf+giant stars, PMS 
stars, and Class I candidates, respectively. The dwarf and giant locus were from Bcsscll & Brett (1988), and T Tauri locus was from Meyer, Calvet, & 
Hillenbrand (1997). The reddening vector was calculated using the reddening law from Nishiyama et al. (2006) and transformed from IRSF to 2MASS 
system. Right panel: H — K s — 0.15 vs. Pjj diagram, dashed and dot-dashed lines are the relationship for interstellar molecular clouds and its upper limit ^ 
(Jones 1989). Solid line shows the result of linear fitting to all data points. Dotted line shows the mode value of Ph- ^ 
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Figure 7 — Histogram of Ph/Ah for all the stars except for protostar candidates. Dashed line on the left side shows 
Ph/Ah = 1.69, which corresponds to the solid line in Fig. 6b. Dashed line on the right side shows Ph/Ah = 3.38. 
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Figure 8 — Relationship between stellar mass and Ph /Ah in NGC 2024. The value of stellar mass was taken from 
Levine et al. (2004). The mass range covers from brown dwarfs to solar-type stars. The vertical and horizontal 
broken lines denote the upper limit of brown dwarf mass (0.08 Mq) and Ph /Ah = 3.38, respectively. 
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Figure 9 — Histograms of Ph for the stars with Pjj > 3 % (left panel) and Ph < 3 % (right panel). Dark gray, light 
gray and white colors show the number of protostar candidates, PMS stars, and dwarf+giant stars, respectively. 
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Figure 10 — Polarization vector of each stellar source on H intensity image. Filled gray circles, open diamonds, and 
open circles denote dwarf+giant stars, PMS stars, and protostar candidates, respectively. 
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Figure 11 — Polarization vectors obtained through the near-infrared polarimetry at H (yellow, this work), dust 
emission polarimetry at 850 fim (blue, Matthews, Fiege, & Moriarty-Schicvcn 2002), and at 100 (im (red, Hildcbrand 
et al. 1995; Dotson et al. 2000). Vectors from the dust emission polarimetry were rotated by 90 degrees for the 
comparison with the E- vector of dichroic polarization (inferred direction of magnetic field). The background image is 
the 850 ^m dust continuum intensity map (logarithmic scale) kindly provided by Doug Johnstone. The contour of 
gray broken line denote the level of 0.9 Jy beam -1 . Plus symbols denote far-infrared source, FIR 1-7 from north to 
south, reported by Mezger et al. (1988, 1992). 



